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Abstract 

A receding-front model for drying of porous material is proposed that explains their drying-rate 
curves based on the dynamics of the evaporation front. The falling-rate regime is attributed to 
the slowing down of the front's propagation inside the medium due to the resistance offered by 
the disorder generated by porosity. The model is solved numerically and the resulting drying-rate 
curve is obtained for the falling-rate period. The curve shows a linear behavior at early times in 
conformance with experiment. 
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I. INTRODUCTION 



Drying refers to the process of removal of free moisture from the surface as well as the 
interior of a solid material by the application of heat. The output from a drying experiment 
can be plotted as a 'drying-rate curve', which is a graph of the drying rate, Q = —A^dw/dt 
(moisture loss per unit time, t, per unit area, A, of the exposed drying surface), versus 
moisture content, w (moisture mass per unit mass of the dry solid) of the material. For 
porous material in general, a common feature of such curves is the existence of two well- 
defined regimes, namely, a constant drying-rate period succeeded by a falling-rate stage that 
is at first linear . The constant-rate period, wherein a film of free moisture is always 
available at the drying surface, is governed fully by the transport conditions external to the 
material. At the so-called critical moisture content that marks the end of the first regime, the 
drying rate begins to fall since the interior moisture, due to internal transport limitations, 
can no longer migrate at sufficient rate to the drying surface to keep it saturated. The 
evaporation front, thence, begins to recede from the drying surface into the solid; moisture 
reaching the surface via migration through the pore volume. Thus, during the falling-rate 
period, the internal transport mechanisms govern the drying process. 

Numerous contributions to the theories of drying rate, with particular emphasis on the 
falling-rate regime, can be found in the literature (see □ ODD and the references therein, 
in particular [3[ which provides an extensive bibliography of references dealing with drying 
models for porous solids). The present work introduces yet another drying model that lies 
in the category of the so called receding- front models (see also 0]). However, rather 
than considering the heat and mass transfer as is customarily done, we shall take a different 
physical approach based on the dynamics of the evaporation front according to the following 
picture. The evaporation front defines a moving interface between the wet and dry regions. 
It is formed by the moisture film surface all through the first stage of drying. Consequently, 
as long as the external conditions remain unchanged, the velocity of the interface and hence 
the drying rate will be constant. However, as the front penetrates inside the porous medium 
during the second stage, the disorder in the medium acts as an inhomogeneous friction force, 
opposing the overall advance of the interface. The slowing down of the evaporation front 
propagation, thence, results in the fall of the drying rate. The disorder, generated by the 
randomness of the pore space and the solid matrix within the porous material, does not 
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change with time and, therefore, can be represented as a quenched noise. Quenched noise 
has been widely applied to model such disorder and its nontrivial effect on the motion and 
morphology of an interface is well known j?| 0, 0, ^ , U | . 



After introducing the drying model in the next section, we solve it numerically to obtain 
the drying-rate curve for the falling-rate stage. The curve shows a linear behavior at early 
times in conformance with experiment. The linear falling-rate regime is described in terms 
of the two free parameters of the model that mimic the independent effects of the porosity 
and the external conditions of the material. 



II. THE DRYING MODEL 

As long as a moisture film is available at the drying surface, the evaporation front will lie 
at its free interface. Consequently, provided external conditions do not change, the velocity 
of the interface and hence the drying rate will remain constant. This trivially reconciles with 
the constant-rate period of the drying-rate curve. Let us, therefore, focus on the falling-rate 
period wherein the front forms a rough interface as it penetrates through the medium 

We consider a preparation for the drying material, wherein the exposed surface is flat 
and the cross-sections parallel to this surface are uniform (FIG. [TJ. Take the plane of the 
exposed surface as the xy-plenae and let z = h(x., t) denote the height of the evaporation front 
below the surface at position x = (x, y) and time t. Consider an areal element of the front 
whose projection onto the xy-plane is d 2 x = dxdy. As the front recedes inside the medium, 
this element sweeps out a volume per unit time given by (d t h)d 2 x. This corresponds to a 
pore volume per unit time equal to e(x, z = h)dthd 2 x, where, < e(x, z) < 1 represents the 
local porosity (pore volume per unit total volume) of the medium. Of course, due to random 
fluctuations, the evolution of the front is stochastic; that is, each repetition of the experiment 
from the same initial conditions may result in a different final configuration for the interface. 
Let us, therefore, denote the average over the ensemble of all such configurations by angular 
brackets. Since moisture resides inside pore volume, the rate of change of the moisture 
content will, thus, be given by 

w = — -jrr / < e(x, h)d t h > d 2 x, 
M Jr 

where, p is the density of moisture (water), M is the mass of the dry solid and the integral 
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is over the whole exposed region, R. Hence, we have for the drying rate, that, 

Q=-?- [ < e (x, h)d t h > d 2 x. (1) 
MA Jr. v ; w 

As a first approximation in treating the disorder, which is generated by the (random) spatial 
variation of the porosity, we shall take the inhomogeneity into account only in so far as it 
affects the dynamics of the front (i.e., d t h). That is, we write (JIJ as 

« = i < 2 > 

where, e is now the mean porosity and 

v(t) = — [ <d t h> d 2 x, 
A Jr 

is the average velocity of the interface. Equation (J2J) is the central equation of our model, 
which relates the drying rate to the propagation velocity of the evaporation front. 

To find the propagation velocity, we need to model the time variation of the interface 
height at any position x. The propagation of the front is induced by the flux of moisture 
particles escaping the interface from random positions. This random flux causes the height 
to increase locally at a rate that fluctuates around some average value V. Because of the 
uncorrelated nature of the randomness in the escape process, these fluctuations are described 
by the 'thermal noise', r) t (x.,t), which is defined through 

<7fc(x,t)>=0, <r h {x,t)r ]t {x.',t')>=2D t 8 2 {x-x.')5{t-t'). (3) 

On top of this, we must add the effect of the disorder generated by the random inhomogeneity 
of the porosity of the medium, which directly affects the local velocity of the interface by 
acting as an inhomogeneous friction force. Hence, the dynamics can be described by the 
equation 

= V + th(x, t) + % (x, z = h), (4) 

where, the quenched noise, r/ 9 (x, z), represents the disorder as mentioned in the introduction 
and is defined through 

< 77, (x, z) > = 0, < 77,(x, z)77,(x', z') > = 2D q 5 2 (x - x.')5(z - z). (5) 

Of course, in (JSJ), the averages (not to be confused with < ... > and hence the notation) are 
over the different realizations of the randomness; for a given material, the quenched disorder 



is to be represented through fixing a particular realization. Whence, as far as the quenched 
randomness is concerned, the evolution of the interface is deterministic, i.e., restarting from 
the same initial conditions always results in the same final configuration. The stochastic 
element enters the dynamics solely through the thermal noise. 

Equation (J3J) is by far the simplest equation describing the motion of the evaporation 
front in a porous medium and together with (J2J), therefore, defines our drying model for 
the (linear) falling-rate regime. These two equations yield the following for the drop in the 
drying rate: 

AQ(t) = Q(0) - Q(t) = l R {rj q (^ 0)- < ^(x, h(x, t)) >}d 2 x, (6) 

where, t — represents the start of the falling-rate regime (corresponding to a flat interface 
at the drying surface, i.e., /i(x, 0) = 0). Also, because < d t h >= d t < h > from equation (6), 
(J2J) can be integrated at once to obtain the moisture loss according to 

Aw(t)=w c -w(t) = ^AH, (7) 

where, w c is the critical moisture content and 

Hit) = 4 / < h(x,t) > d 2 x, 
A Jr 

is the average height of the interface. Thus, at the start of the falling-rate regime w has the 
initial value w c , while decreasing as the time elapses at a rate which is directly determined 
by the velocity of the evaporation front. Once the configuration of the interface at arbitrary 
time, t, is determined from its equation of motion, (jUJ) and (J7|) can be used to obtain the 
drying-rate curve as a graph of AQ versus Aw. 



III. NUMERICAL SOLUTION AND RESULTS 



We introduce the dimensionless variables 

1 _ * ~_ * 

X= 4 X ' Z = V t = W 

where, the characteristic length and time scales are defined by the two coefficients entering 
the equation of motion according to the combinations 



These two independent scales are set by the disorder in the medium (through D q ) as well 
as the external conditions that influence the escape process (through D t ); most notably 
the temperature of the hot environment. They remain free parameters of the model to 
be determined ultimately from comparison of the results with experiment. For numerical 
purposes, it is useful to introduce the dimensionless quantities 

1 ~ V ~ 1 1 

M x > *) = T h (*> *)' V = 77' ^( x ' *) = 77^( x ' *)' ^( x ' z "> = T7^( x ' z )' 

tO Vq Vq Vq 

where, 

Vo = T - A' 

is the characteristic velocity. Using Q and ([5]). the second moments of the noise variables, 
therefore, translate into 

< ^(x,f)^(x',f) >= 2 5 2 (±-±')5(t-t'), < fj q (%z)fj g (±',z') > = 2 5 2 (±-±')5(z-~z'), (8) 

while the first moments remain zero, of course. The stochastic equation of motion, thus, 
becomes 

®^ = V + fkfr%+U^)> (9) 
with h(5t, 0) = 0, as the initial condition. To integrate this equation numerically, we dis- 
cretize it by taking the points (x, z) to belong to a cubic lattice of size L 3 and use the Euler 
algorithm with time steps At. (The lattice points should be much more closely spaced in the 
z direction than in the other two directions. This is because the z coordinate of the quenched 
noise in (jUJ) is pre-determined by the value of h.) The noises fj t (% t) and fj q (5c,z) are both 
chosen from a uniform distribution of random numbers between —1 and +1; a distribution 
that, clearly, satisfies the above conditions on the moments. However, the quenched noise 
is generated only once and for all (corresponding to fixing a particular realization of the 
disorder, as pointed out before), while the thermal noise is regenerated at every repetition 
of the calculation. Hence, we generally end up with a different value for /i(x, t) at each 
repetition. 

In terms of the rescaled (dimensionless) rate drop and the rescaled moisture loss, 

AQ(t) = ^-AQ(t) 7 Aw(t) = ^-Aw(t), (10) 
peVo pe£ A 
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respectively, equations (0) and (J7J), thus, reduce to 

AQ(t) = L- 2 ^{r/ 9 (5c,0)- < fj q (x,h(%t)) >} 

X 

Aw(t) =L- 2 Y^< h(i, t) > (11) 

X 

In FIG. |2J we plot the rescaled drying-rate curve as a graph of AQ versus Aw for various 
values of V and L = 40. (Variation of L did not produce noticeable change.) It is seen 
that the curves all collapse into a single straight line at early times, thereby establishing a 
'universal' linear behavior. Using (fTUj) one can, therefore, fit any experimental drying-rate 
curve during the linear falling-rate period by appropriate adjustment of the independent 
parameters £q and to- Our model, thus, describes the linear falling-rate regime in terms of 
two free parameters that represent the independent effects of the porosity and the external 
conditions of the material. 
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FIG. 1: Evaporation front in porous medium 
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FIG. 2: Rescaled drying-rate curve in the falling-rate regime 




